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Abstract. Recent observations of two black hole candidates (GX 339-4 and J1753. 5-0127) in the low-hard state (Lx/i-Edd - 
0.003 - 0.05) suggest the presence of a cool accretion disk very close to the innermost stable orbit of the black hole. This runs 
counter to models of the low-hard state in which the cool disk is truncated at a much larger radius. We study the interaction 
between a moderately truncated disk and a hot inner flow. Ion-bombardment heats the surface of the disk in the overlap region 
between a two-temperature advection-dominated accretion flow and a standard accretion disk, producing a hot (kT^ = 70keV) 
layer on the surface of the cool disk. The hard X-ray flux from this layer heats the inner parts of the underlying cool disk, 
producing a soft X-ray excess. Together with interstellar absorption these effects mimic the thermal spectrum from a disk 
extending to the last stable orbit. The results show that soft excesses in the low-hard state are a natural feature of truncated disk 
models. 
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1. Introduction 

The geometry of the low-luminosity ("low-hard") state of 
Galactic Black Hole Candidates (GBHC), in which the spec- 
trum is dominated by a power law X-ray flux extending to high 
energies, has been an open question for several decades. While 
it is generally believed that the power law spectrum is formed 
by inverse Compton scattering, there is no consensus about the 
geometry of the flow, source of seed photons or energy distri- 
bution for the Comptonizing electrons. 

Broadly speaking, there are two classes of model to explain 
the spectrum in the low-hard state. The first is the "corona" 
model, in which the disk remains untruncated or nearly un- 
truncated at luminosities Lx - lO^^'Ledd- The hard power law 
spectrum comes from a hot and patchy corona (perhaps pow- 
ered by magnetic flares (Idi Matted 1 1998t lBeloborodovlll999t 
Merloni & Fabian 200 lb) on top of the disk, while the sur- 
rounding region is bombarded with high energy photons, pro- 
ducing the observed reflection and Fe-K fluorescence compo- 
nents. In the alternate, "truncated disk" model the thin disk is 
truncated at some distance from the black hole and the inner 
region is filled with a hot, radiatively-inefficient flow, which 
produces the hard spectrum. The reflection spectrum and Fe-K 
fluorescence is then produced by the interaction of the hard X- 
rays with the inner part of the truncated disk, or in some cool 
outflow moving away from the disk. For a recent d iscussion of 
the low-hard state see sect. 4 of iDone et al.l OOOTh . 

In theory, the presence or absence of a cool disk should be 
confirmable through direct detection of a soft X-ray blackbody 



component at low energies. In practice however, this is made 
difficult by the fact that at low accretion rates the temperature 
of even an untruncated disk will drop from about l-2keV in the 
high soft state to ~ 0.1-0.3keV, which puts it out of the range 
of most X-ray detectors. Additionally, the effects of interstellar 
absorption become very strong at around 0. IkeV, so that detect- 
ing a soft excess and accurately measuring its parameters will 
depend somewhat on how accurately the interstellar absorption 
can be determined. 

Even with these challenges, a soft excess in the low-hard 
state has previ ously been reported in several sources. The first 



was Cyg X-1 jBalucinska-Church et all 1 19951: iDi Salvo et al 



20011) . although its association with an accretion disk is com- 
plicated by the fact that Cyg X-1 is a high mass X-ray bi- 
nary accreting from a wind. This question was also the focus 
of two recent papers, iMiller et al.l (l2006allbl) . in which the au- 



thors studied long-exposure XMM-Newton spectra of two dif- 
ferent GBHCs, SWIFT J1753.5-0127 and GX 339-4, at low 
luminosities (Lx/^Edd ~ 0.003 - 0.05). Soft excesses at simi- 
lar luminosities in these tw o sources have also b een reported in 



Ramadevi & Seethal ( l2007h (J1753.5-0127) and lTomsicketal. 



(120081) (GX 339-4). Since these two observations, there have 
also been observations of soft excesses in several other sources. 
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Rvkoff et al.l (I2007h made several observations of the soft com- 
ponent of XTE J1817-330 with Swift during the outburst de- 
cline of that source down to a luminosity of Lx/Lem ~ 0.001, 
whil e a soft component in G RO J 1655-40 has been repor ted by 
both Brocksop p et al.l (120061) and lTakahashi et al.l (l2008h using 
different telescopes. 
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To inte rpret the soft excesses i n SWIFT J1753. 5-0127 and 
GX 339-4, iMiller et al.1 (l2006albb fit the data with a several 
XSPEC models, trying various black-body disk models and 
simple hard X-ray components (both a power law and various 
Comptonization models). In GX 339-4 a broad Fe-K line was 
also observed and fit with a relativistically broadened reflec- 
tion model. Using blackbody models for a standard accreting 
disk, the authors found disks with maximum temperatures of 
kT ~ 0.2-0.4 keV, and inner radii consistent with the innermost 
stable circular orbit of a black hole. 

At the inferred low accretion rates in the hard state, a disk 
extending to the last stable orbit would produce a soft X-ray 
component with peak close to the cutoff due to interstellar ab- 
sorption. Unless an accurate independent measure of the in- 
terstellar absorption column is available, spectral fitting pro- 
cedures cannot reliably distinguish between a thermal peak at 
kT - 0.3 keV with one interstellar absorption column and a 
cooler component with a lower energy component cutoff by a 
slightly higher interstellar absorption column. 

For energetic reasons the hard X-ray component which 
dominates the luminosity in the hard state must originate near 
the black hole, the same region as the proposed cool disk. Some 
form of interaction of hard X-rays with the cool disk must take 
place, and this implies that the isolated cool disk models used 
as 'components' in fits to observed spectra are unrealistic. In 
fact, most models for the hard X-ray component include some 
prescription for the reprocessing of hard into soft radiation, 
whether these be truncated disks or extended disk models. As 
shown by Haardt and Maraschi (1991), such models generi- 
cally produce a similar energy flux in soft and hard X-rays. A 
strong soft component is thus a natural consequence in trun- 
cated as well as extended disk models for the hard state. 

The main difference in a truncated disk model is that the 
soft flux originates from a larger surface area and consequently 
has a lower temperature, putting its spectral peak below 0.5 
keV. After interstellar absorption the soft component has a peak 
around 0.5 keV that can be mistaken for an apparent thermal 
peak with the temperature of a disk near the last stable orbit. 

In this paper we examine this with a more quantitative 
model for truncated disks. At the inner edge of a truncated disk 
the accretion flow must change in nature from a relatively cool, 
thin disk into a much hotter, vertically-extended inner flow. 
There will thus necessarily be so me interaction between th e 
two, either through radiation (e.g . iHaardt & Maraschil (1199 11) ) 
or matter exchange (ISpruitil 1997b . or both. Our goal is to deter- 
mine whether such a mo del could reproduce the soft spectral 
components reported by iMiller et al.l (l2006allbl) . We will find 
that disks truncated at 15-20 Schwarzschild radii can in fact 
produce soft components of the observed strength and shape. 
An alternative model i nvestigating re-co ndensation fron am 
ADAF is considered bv lTaam et al.l (l2008b . 



of how luminosity in the disk is produced. In the higher lumi- 
nosity "high-soft" state, (approximately) blackbody flux from 
the disk dominates the spectrum. The disk's inner radius can 
be inferred, provided the distance and inclination of the source 
are known, from the assumption that the radiation is produced 
by internal viscous dissipation in the disk as it accretes, and an 
assumed 'color correction' to the blackbody spectrum. 

However, when the spectrum becomes dominated by very 
hard X-ray and y-ray radiation, the radial temperature and lu- 
minosity profile of the disk will depend on the interaction be- 
tween the hard radiation and the disk. This is because most of 
the accretion power is now in the hot gas producing the hard 
radiation, some of which will interact with the cool disk (as 
is seen in the reflection spectrum and Fe fluorescence), and 
may provide a substantial source of heating. This heating of 
the disk surface converts the hard radiation into a soft compo- 
nent. I n the simplest version of this model, Haardt & Maraschi 
(Il99lh calculated the energy balance for a hot corona covering 
a cool disk and found that the flux in the soft and hard com- 
ponents will be roughly equal. This predicts that a substantial 
soft X-ray component is a universal feature of the hard-state 
spectrum, whether the disk is truncated or not. Its detectability 
depends on the sensitivity of detectors in the 0.1-1 keV range, 
and the interstellar absorption column. 

For models in which a cool, truncated accretion disk en- 
circles a very hot inner flow (which cools through inverse 
Compton scattering), there will necessarily be a region of inter- 
action near the truncation radius, which will heat the inner edge 
of the disk. This will come either from hard photons bombard- 
ing the disk as in the Haardt-Maraschi model, or through matter 
interaction, with hot protons from the inner flow directly collid- 
ing with the cool disk. In both cases previous work has shown 
that a soft component is produced, but a more detailed model is 
needed to predict the radiation spectrum. The first case, exam- 
ining the structure of a thin accretion disk bo mbarded by hard 



photons ha s been studied by for AGN di sks iBallantyne et al 



2. Physics of Interaction Region 

2.1. Origin of Soft Excesses 

Determining the inner radius of an accretion disk from its spec- 
trum relies on (among other things) an accurate understanding 



(120011) and Nayakshin & Kallman (2001), and more recently 
for GBHC disks bv .Ross & Fabian (.2007) . In this paper we fo- 
cus on the second case, in which a moderately truncated disk is 
embedded in a hot, two-temperature advection-dominated ac- 
cretion flow (ADAF). 

2.2. Definition of tlie model 

We begin with the results from a prototype model for ion 
bombardment on cool disks, initially propos e d by Spruii 
( 1997), and extended in Spruit & Deufej (l2002l) . iDeufeletal 
(2002) and lDullemond& Spruit (2005r (hereafter DS05). In 
this model a cool disk is embedded in a two-temperature ADAF 
in which the protons are close to their virial temperature (~ 20 
MeV). The protons bombard the disk, and are stopped via 
Coulomb collisions with the disk's electrons, thus transferring 
their energy to the disk. The energy from the protons is suffi- 
cient to evaporate the upper layers of the disk into a hot corona 
with kTc -60-80 keV (called the "hot layer" in our nomen- 
clature, although the temperature in this layer is still much 
cooler than the virial temperature), whose temperature is set by 
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Fig. 1. Schematic structure of a cool disk embedded in an 
ADAF. Here the ADAF (the transparent outer region) extends 
over the inner edge of the cool Shakura-Sunyaev disk (SSD), 
shown by the solid thin disk, bombarding it with high energy 
protons {kTp ~ 20 MeV). This evaporates the surface of the 
disk into a (kT^ ~ 70keV) corona, called here the hot layer 
(HOT LAYER, hatched). The higher viscosity of the hot layer 
causes it to spill over the inner edge of the SSD where cool- 
ing through inverse Compton scattering is no longer efficient 
(HOT RING), and thermal instability causes the layer to evap- 
orate into the ADAF. The mass transfer is represented by the 
small arrows. 



a balance between heating from the ions and cooling, predom- 
inantly through i nverse Compton scattering of disk photons. 
Deufel & Spruit! (120001) found that the optical depth of the hot 
layer is around unity, varying only weakly with irradiating flux 
and distance from the hole. 

The higher viscosity of the hot layer also allows it to 
spill over inside the in ner edge of the coo l disk (see fig. [T]i- 
Using 1-D simulations. iDeufel et al.l (120021) suggested that the 
lack of seed photons in this layer will cause the tempera- 
ture to rise to about kTe ~ 200 - 300 keV, since cooling 
from inverse Compton scattering will be much less efficient. 
Spruit & Deufell (120021) found that this region will become un- 
stable and evaporate into the ADAF. This shows how an ADAF 
can be maintained inside a truncated cool disk. The key for the 
whole process is the presence of a component of intermedi- 
ate temperature (the hot layer). On one hand this component 
produces a hard, Comptonized spectrum, while on the other its 
evaporation feeds the ADAF. 

DS05 extended this work to examine the radial depen- 
dence of the three components of the flow. They set up a 
model in which the truncation radius of the cool disk is a 
free parameter, and the three layers are assumed to interact 
only via the mass transfer processes outlined above. They as- 
sumed the cool disk was a standard Shakura-Sunyaev disk 
dShakura & Sunvaeviri973l) . and that the ADAF flo w temper- 
ature followed the model of iNaravan & Yil (119941) . They ap- 
proximated the temperature in the hot layer to be constant, and 
also assumed a constant temperature for its extension inside the 
inner edge of the cool disk (which we will call here the "hot 
ring"lj. With these assumptions, the radial profiles of the three 



components can then be determined by solving the equations 
for mass and angular momentum conservation for a thin disk. 

This simplified analysis yields for each component esti- 
mates of the mass accretion rate M,, heating rate Qi (from in- 
ternal viscous dissipation and ion heating), and surface density 
E, as a function of radius. For further discussion of the model, 
see sect. 2 of DS05o The input parameters for the model are 
the black hole mass M^ff, the magn itude of the a viscosity pa- 
rameter dShakura & Sunvaevil 1973b . and the truncation radius 
RssD for the disk. A sketch of the model with its various com- 
ponents is seen in fig. [T] Assuming we know the dominant ra- 
diation mechanisms, we can determine the spectrum from the 
resulting energy distribution for each layer. This is the primary 
goal of the current work. 

2.3. Energy and mass balance 

The three components of the model exchange mass and energy. 
In the steady state assumed here, the sources and sinks bal- 
ance in each. The topmost layer, the ADAF, loses mass and 
energy to the hot layer, at rates per unit area of the order Cs;,pi 
and CsiWi, respectively. Here Cj^, pi, and wi are respectively the 
thermal speed in the ADAF (i.e. near the virial speed or orbital 
velocity), the density, and enthalpy. Depending on its radiative 
efficiency the ADAF could make a substantial contribution to 
the hard X-ray flux in addition to the flux from the ion-heated 
hot la yers. For our prot otype model we ignore this contribution, 
since Esin et al.l (Il997h estimate the efficiency of an ADAF to 
be about e ~ ma^ ^ 1.5 x lO"'* (where in is the accretion 
rate scaled in terms of Eddington). If this is the case, the lu- 
minosity of the ADAF will scale with m^, while the energy 
influx into the cool disk scales as p-, ~ m, (since T ~ Tvii in 
the ADAF), meaning that the ADAF's indirect contribution to 
radiation through heating the cool disk will dominate at low lu- 
minosities. However, the radiative efficiency of hot inner ring 
is strongly model dependent, and may be much higher, which 
would make the spectral contribution of the ADAF important. 
We discuss this question further in sect. H] 

There is a mass flux by evaporation from the hot layer into 
the ADAF, as derived in I Spruit & Deufell (120021) . The energy 
flux into the ADAF associated with this evaporation can be ne- 
glected, since the temperature of the hot layer is low compared 
with the temperature of the ion flow. The hot layer is heated in 
similar measure by viscous dissipation and the hot ions it ab- 
sorbs from the ion flow above (Qrac.w = 2 x 10'^ ergs^'cm"^, 
and 2,;w = 7 X 10'^ ergs^'cm"^, respectively for the reference 
model). It is cooled by Compton upscattering of disk photons 
passing through it, which will produce a hard X-ray spectrum. 

There is also a mass flux feeding mass into the hot layer 
from the cool disk underneath; it is parametrized as in sect. 
2.2 of DS05. The cool disk is heated in three ways: by internal 
viscous dissipation, by reprocessing of hard (Comptonized) X- 



' The nomenclature we adopt in this paper is slightly different from 
DS05. They referred to the ADAF as an "ISAF", or Ion-Supported 
Accretion Flow, while what we have called the "hot inner ring" they 



term the "hot layer", and our "hot surface layer" is in that paper re- 
ferred to as the "warm layer". 

^ We note that in that work the numerical factor of equation (2) is 
incorrect, it should be 2.64 x 10^, and equa tion (3) should be the same 
as equation (16) in lSpruit & Deufell ( l2002h . 
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rays from the hot layer above, and by the (small) number of hot 
ions from the ADAF that pass through the hot layer. 

The reprocessing of the hard X-rays is further divided into a 
fraction a that is reflected, and a fraction (1 -a) that is absorbed 
and assumed to thermalize completely. The reflection process 
and the spectrum resulting from it are based on results from 
the literature; this is described in sect. 12.4.21 The spectrum of 
the cool disk feeding photons to the Comptonization process 
in the hot layer thus consists of a blackbody component plus a 
reflection spectrum. 

The surface temperature of the cool disk, assumed to be 
close to its effective temperture, is given by 



o-bT, 



eff 



Gc,vi,sc + (l-fl)/^h +/cfii- 



(1) 



where 2c,visc is the viscous heating rate of the cool disk c, a its 
X-ray reflection albedo, F^ the flux of downward directed X- 
ray photons from the hot layer h, and f^Qi the fraction /c of the 
energy flux Qi in hot ions that make it through the hot layer and 
are absorbed instead in the cool disk. The albedo a of the disk 
depends on how ionized the disk is, which we discuss more in 
sec. l2A2l 

Compton upscattering of the soft flux from the cool disk 
determines the detailed shape of the spectrum emerging from 
the hot layer. In addition there is an X-ray component from the 
'hot ring' interior to the inner edge of the cool disk. Since it 
receives fewer of the soft photons and hence is hotter, it adds a 
harder component to the spectrum. 

These ingredients are the same as in DS05, except that 
the reflection process is included more realistically, and the 
Comptonization of soft photons is treated in detail with the 
Monte Carlo code of Giannios (see sec. O so a realistic X-ray 
spectrum is obtained. 

2.4. Reference model 

We illustrate the model first with representative parameter val- 
ues for a black hole binary. In sects 13. Il and l3.2l the parameters 
are adjusted for application to observations of specific objects. 

We take Mbh - IOMq for the mass, Ri„ - 20 Rs for the 
truncation radius (where Rs - 2GM/c- is the Schwarzschild 
radius). A typical if somewhat large viscosity parameter a = 
0.2 is used for all accretion components. The flow model of 
DS05 yields an accretion rate of M = 3.1 x 10"' g/s for these 
parameter values. In terms of the Eddington accretion rate (de- 
fined here with an assumed bolometric efficiency of 10%) this 
translates to M = 2.2 x lO^^MEdJj, and bolometric luminos- 
ity is Lx - 4.2 X 10"'*LEdd, which is significantly lower than 
observed luminosities for the low/hard state. We discuss this 
discrepancy and possible solution further in sects [3] and |4] 

In fig. |2]we show the temperature profile of the cool disk 
component of this reference model. The dashed line shows the 
temperature profile that would result if only the viscous dissipa- 
tion in the cool disk were included. The cool disk is completely 




^ In DS05 a plotting error resulted in the mass accretion rates in 
all figures being overstated by a factor of 100 with respect to the 
Eddington rate. 



1000 



Fig. 2. Eff'ect of ion-bombardment on the surface temperature 
of a cool (SSD) disk with inner edge at 20Rs (for the other 
model parameters see text). Black solid line: radial tempera- 
ture profile when disk heating is considered. Green dashed line: 
radial temperature profile in disk without considering heating 
from hot layer Blue dotted line: radial temperature profile for 
an untruncated disk at the same accretion rate. 

evaporated into the hot layer before the inner edge of the accre- 
tion disk, hence its surface density drops to zero at the inner 
edge of the disk in the same way as for a standard cool disk ac- 
creting on a slowly rotating object (the standard "zero-torque 
inner boundary condition"). If there were viscous coupling be- 
tween the inner edge of the disk and the hot ring, however, the 
resulting temperatures near the inner edge would be higher than 
in our current model. This would result in a stronger soft com- 
ponent, whose contribution would remain observable for larger 
truncation radii than we have considered here. The solid line 
shows the efi'ect when heating from the hot layer is included. 
As can be seen in the figure, the eff'ect of heating extends to 
large radii. The dotted line shows the temperature profile for an 
standard untruncated disk for the same accretion rate. 

We can also compute the bolometric radiative efficiency 
for the model, assuming that the hot layer and inner ring ra- 
diate efficiently. For our reference model with a disk truncated 
at 20 Rs, we find a radiative efficiency (the ratio of the total 
luminosity to the accreting rest-mass energy flux) of about 4%. 
In comparison, the radiative efficiency of a disk truncated at 
20 Rs alone is about 1%, which demonstrates that the ion- 
bombardment process can significantly increase the radiative 
efficiency. 

2.4.1 . Cool disk emission 

We begin by estimating the spectrum of the radiation produced 
by the thin disk, i.e. the seed photons of the Comptonization 
in the hot layer. At large effective optical depths in the disk 
the radiation field is a blackbody. The flux emerging from the 
photosphere of the disk is reduced, however, by Thomson scat- 
tering, resulting in a 'modified blackbody', a spectrum with 
same temperature but reduced flux. In the literature, it is of- 
ten described equivalently by a 'color correction' factor /coi. 
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Fig. 3. Surface density (2) profile for reference model, where 
Mbh = lOMo, a = 0.2 and R,,,^ = 20/?,. Green dashed 
line: Shakura-Sunyaev disk. Blue dotted line: Hot layer. Orange 
dash-dotted line: Hot ring. Red solid line: ADAF. 



measuring the effect instead as the ratio of color temperature 
to effective temperature, Tco\ - /coiT^eff. Several papers have 
done extensive calculations of this spectral hardening in the 
outer layers of a disk dom inated b y electron sca ttering, in- 
cludingls himura & Takahara ( 1995), Merloni et all fcoOO). and 
Davis et alj (l2005b . and concluded that the the color-correction 
is of the order f^d - 1 .4 - 2, with the magnitude of correction 
increasing gradually as the luminosity increases. The approxi- 
mation of a single parameter across the entire spectrum become 
less good as the flux decreases. 

In a cool disk the (vertically integrated) viscous dissipation 
Qvisc is balanced by radiative loss at its surface, 



2visc - O-BT^fi - O- b(T col /f col) ■ 



(2) 



The value for f^oi will change the determination of the inner 
disk radius, since the luminosity will be depressed by a factor 
of f^J^ relative to the temperature, so that when color correc- 
tion is considered, the inner radius of the accretion disk will be 
increased by a factor f^^y 

In this work we have the additional complication of the 
hot surface layer, which will heat the upper layers of the 
disk substantia l ly an d may change the spectral hardening. 
Ross & FabianI (l2007h have estimated the effects of incident 
hard X-ray radiation on galactic black hole disks, although they 
focus on systems in which the disk is hotter (kT ~ 0.35 keV) 
than ours and dominates the overall spectrum. However, they 
do find substantial spectral hardening, and disk temperatures 
which increase with increasing X-ray flux (fig. 4 of that paper). 
They also find that the incident radiation can change the ion- 
ization structure of the upper layers of the cool disk, so that the 
shape of the disk spectrum is significantly altered by absorption 
lines (see their fig. 6). 

Since it is beyond the scope of this work to do a similar 
radiation transfer calculation, we instead choose a single color 
correction factor, fc„i = 1 .7 as a typical value for all our mod- 
els. 



2.4.2. Reflection spectrum 

The hard flux incident on the cool disk undergoes reprocess- 
ing by electron scattering and atomic processes, giving rise to a 
series of emission lines (especially the Fe-K fluorescence line) 
and a Compton reflection hump. In the course of these pro- 
cesses, part of the incident energy flux ends up heating the 
plasma. Detailed calculations of these processes is beyond the 
scope here. Instead, we simplify the physics by separating the 
incident flux into a reflected part and a part that is treated as be- 
ing absorbed. We treat the absorbed part in the same way as the 
energy input by viscous dissipation in the disk, i.e. a fraction 
(1 - a) is added to the left hand side in (|2]i. 

The remaining part is treated by a modi fied version of 
the X SPEC model "REFLION" described in JRoss & Fabian 
2005b . This model calculates the reflection spectrum of a power 
law flux between 0.1-300 keV hitting a constant density AGN 
disk. We set the lower cutoff of this power law at 0.2 keV, since 
we expect disk temperatures around 0.1 - 0.2 keV. REFLION 
calculates the reflection spectrum for a fixed ionization param- 
eter, ^ = AuFxri^ , the ratio of the incident flux on the disk to 
the hydrogen number density. The disk, however, is stratified, 
so an estimate has to be made of the depth in the disk atmo- 
sphere where most of the reprocessing takes place. We assume 
here that this depth is just the electron scattering photosphere 
of the cool disk, Tes = 1 . By hydrostatic balance, the pressure 
at this depth is approximately 



P^ — , 



(3) 



where g is the vertical component of the acceleration of 
gravity at the height of Tes - 1 above the midplane. Assuming 
this height is about twice the nominal disk thickness H = 
Cs/Qk, the pressure is P ^ 2csQk- Here Cj is the sound speed 
at the midplane of the disk, which is determined from the heat- 
ing rate gvisc by ^ standard thin disk model. To fix the density 
corresponding to this pressure the temperature has to be deter- 
mined. We assume for this the effective surface temperature T^ff 
that follows from the viscous heating rate. If F^ is the incident 
X-ray flux from the hot layer, the ionization parameter is then 



f- 



0.kCs 



(4) 



The value of ^ varies with distance r from the black hole; 
an energy- weighted mean is ^ ^ 10 erg cm s"' for the refer- 
ence model. We use this as a representative value for the re- 
flection calculation instead of an integration over the disk. For 
this value of ^, the resulting albedo is low (^ 0.2), so that most 
of the hard flux is absorbed and thermalized in the disk. The 
total reflection and disk spectrum (before being passed through 
the hot layer), is seen in the red long-dashed line in fig.m This 
is the input spectrum for the Comptonization by the hot layer. 
Since the optical depth of the hot layer is only around unity, 
it will also make a significant direct contribution to the output 
spectrum. 
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2.4.3. Uncertainties in tine reflection spectrum 

Through its effect on the gas density, the temperature assumed 
in the reprocessing region has a direct effect on the ionization 
parameter and the resulting reflection spectrum. Unfortunately, 
there is a significant uncertainty in this temperature, since it de- 
pends itself on the ionization parameter At low densities higher 
in the atmosphere, the ionization parameter is high and the tem- 
perature will be close to the Compton temperature of the inci- 
dent X-ray spectrum. For our reference model, this would be 
around 3 keV. Deeper in the atmosphere, the ionization param- 
eter would be lower and the temperature closer to equilibrium 
with the energy density of the radiation field, of the order of the 
effective temperature of the cool disk. Since the transition be- 
tween these regimes takes place around the reprocessing layer 
itself, all depends on details of the radiation physics. 



Several groups dBallantvne et al.l 1200 it iNavakshin et al 



2000 !) have calculated the change in the vertical structure of 
AGN disks as a result of incident hard radiation, and found 
the upper layers become stratified, with the top being heated to 
close to the Compton temperature and a sharp transition to the 
inner layer which radiates at the disk's effective temperature. 
They find that this stratification wi ll change the disk's refl ection 
properties substantially. Although lRoss & FabianI (l2007b study 
a similar situation in galactic black holes, and find no simi- 
lar stratification, they focus on cases in which the disk domi- 
nates the luminosity. For lower disk ten iperatures stra tification 
could presumably still occur. However. INavakshin & Kallman] 
(1200 ll) studied the case where a disk was overlaid with a sur- 
face corona, and included the effect of its weight on the gas 
pressure. They found that this was sufficient to prevent stratifi- 
cation. 

An additional caveat in using the constant density mod- 
els of IRoss & Fabian (2005) is that they consider AGN disks, 
which are much cooler and less dense than those in GBHCs. 
As a result, the reflection spectra do not consider the flux in 
the disk itself (which can change the ionization state signifi- 
cantly), and the effects from things like three-body recombi- 
nation, which will become im portant as th e disk densities be- 
come higher Recent work bv IRoss & Fa bian (2007) has ex- 
amined reflection spectra in GBHCs, although they focus on 
systems in which the disk dominates the spectrum and is hotter 
(kTe - 0.35 keV) than the disks considered in this work. 

2.4.4. Comptonization 

To calculate the spectrum from inverse Compton scattering 
through the hot layers of the disk we use a one-dimensional 
Monte Carlo simulation. The radial variation of quantities in 
the disk and hot layer is replaced by representative values, 
for which we use energy-weighted means. The simulation as- 
sumes a slab geometry, with the cool optically thick disk below 
a much hotter surface layer with moderate optical depth. For 
the seed photons, we use the disk+reflection spectrum found 
in the previous two sections. To capture the emission lines 
in this spectrum, we set the resolution of the simulation to 
AE/E = 0.046. The number of seed photons followed through 
the hot layer is of the order 10^, sufficient to represent the re- 



sult to a noise level comparable with typical observations. The 
output spectrum is angle dependent (because of the increas- 
ing optical depth with inclination), so a value has to be as- 
sumed for the inclination angle 6 to the line of sight (measured 
from the normal of the disk). For our reference model we use 
yu - cos 6 = 0.5 as a representative value. 

The input parameters of the Comptonization calculation are 
the input spectrum, the optical depth and the temperature of 
the hot layer. The resulting total energy flux in Comptonized 
photons (integrated over the spectrum and summed over both 
sides of the hot layer) has to match the energy input into the 
hot layer: the sum of ion heating and viscous dissipation. We 
can estimate the optical depth for the hot layer from the surface 
density of the hot layer, where (t) - Kesi^Hot)- The temperature 
of the hot layer is adjusted iteratively in the calculations to meet 
this condition to an accuracy of 10%. For the reference model, 
we find kTc - 70keV, and t = 0.87. The resulting emergent 
spectrum of the hot layer is shown in the blue dotted line in fig. 
m for an inclination angle of /i = 0.5. The photon index of this 
spectrum in the range 2-10 keV is about F = 1 .96. 



2.4.5. Treatment of the hot ring 

In the hot ring (where the hot layer has spilled over inside in- 
ner edge of the cool disk) there is no underlying cool disk any 
more and it receives its soft photons only by scattering from 
larger distances: it is 'photon starved' compared with the hot 
layer itself. Since the energy input by ion heating and vis- 
cous dissipation are still similar, the temperature is higher and 
Comptonization correspondingly stronger The hot ring there- 
fore makes its contribution mostly at the high energy end of the 
spectrum, and is less important for the 'soft component' of the 
spectrum that is the focus of our study. We include it, however, 
since we also want to achieve a reasonable fit to the overall 
spectral energy distribution in the observations discussed in the 
next sections. 

The calculation of DS05 relied on the earlier one- 



dimensional work of iDeufel et al.l (l2002h in order to set the 
temperature and energetic contribution from the hot inner ring, 
which predicts a very small contribution from the hot ring. With 
a sufficiently detailed geometrical model for the hot ring, the 
soft photon input by scattering could in principle be modelled 
more accurately, but the level of detail needed is probably be- 
yond the limits of the present model. We therefore treat the 
soft input flux in this component of the model as an adjustable 
unknown when fitting to spectra. For this, we introduce a pa- 
rameter (, which represents the fraction of seed photons from 
the cool disk that cool the hot layer, so that 1 -( goes to cool the 
hot ring. For the reference model, ( is set by the contribution 
predicted from the DS05 model for the hot inner ring. 

Also treated as adjustable is the temperature reached by 
equilibrium between heating and Comptonization. We assume 
that cooling in this region is still moderately efficient, choosing 
temperatures in the range kT^. = 180 - 200 keV. For simplicity 
(since the angular distribution of seed photons is also uncertain) 
we model the hot ring with a plane-parallel Monte Carlo simu- 
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lation with the same resolution as in the hot layer. We discuss 
the limitations of this approach in sect. |4] 

As in the hot layer, the X-ray energy flux produced by 
Comptonization has to match the energy input by viscous dis- 
sipation and ion heating. Together with the now assumed value 
for the temperature, this determines ^, and the optical depth of 
the ring. For the reference parameters of this section, we find 
^ - 0.99 and an optical depth of ~ 0.7. At a temperature of 200 
ke V we find that the ring contributes only 11 % of the overall 
hard flux (Fx - 2 - 200 keV) for the reference values a = 0.2 
of the viscosity and Ri„ = 2QRs . 

This is because the radial width of the hot ring has a rather 
limited extent, since the evaporation process is very efficient. 
The hot layer evaporates very quickly after it has flowed over 
inner edge of the cool disk. The hot Comptonized component 
from the hot ring is shown in the orange dash-dotted line in fig. 

El 

2.5. Results for the reference model 

The final spectrum and its various components is shown in the 
top panel of fig |4] The accretion rate for the reference model 
is M/Mem - 0.002 assuming an efficiency of 10%, while the 
luminosity in the 0.5-10 keV range is Lx/Z-Edd = 10""*. The 
total spectrum is shown in black. The individual components 
run as follows. The green thermal component (dash-double dot- 
ted line) shows the spectral contribution from the outer part of 
the disk where the hot layer is no longer significant (outside 
R/Rs = 100), while the red long-dashed component shows 
the rest of the disk and reflection spectrum. For the reference 
model the reflection ionization parameter is small enough that 
the reflection and iron line are not apparent in the final spec- 
trum, although we again stress that we are using a reflection 
model developed for AGN, so in reality the reflection could be 
stronger. The blue dotted line shows the Comptonized spectrum 
from the hot layer, while the orange dash-dotted line shows 
the Comptonized spectrum from the hot ring. Fitting the over- 
all spectrum with a photon index F = 1.91 in the 1-10 keV 
range, we see a small soft excess below 0.5keV, even though 
the maximum disk temperature is only O.OSkeV. The bottom 
panel of fig |4] shows the total spectrum divided by a power law 
with F = 1.91. The observed deviation from a single power 
law in the hard part of the spectrum (which leads to a deficit 
around 1 keV and a harder power law index above 10 keV) is 
caused by anisotropic Comptonization res ulting from co nsid- 
ering a plane-parallel configuration (see e.g. lHaardti (I1993l) and 
sect. Hlof this paper), and also by the contribution from the hot 
ring. Except for the very low accretion rat e, we see a spectrum 
that is qualitatively similar to those of .Miller et al.l (l2006a) and 
MmeretalJ(l2006bl) . 



3. Comparison to observations 

The luminosity of our reference model described above, of the 
order lO^^^Lsdd is substantially lower than the luminosities (of 
the order 3 x lO^^-lO^^Lndd) inferred for the observed sources 
to which we want to apply the model. This is a consequence of 
the flow model in DS05 that is the basis of our analysis. In it, 
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Fig. 4. Top: Relative contibution from each component for ref- 
erence model, a = 0.2, Ri„ = lOR^, M = IOMq. Red long- 
dashed line: Spectrum from the disk (modified blackbody plus 
reflection spectrum). Green dash-double dotted line: Spectrum 
from outer disk. Blue dotted line: Comptonized spectrum from 
hot layer. Orange dash-dotted line: Comptonized spectrum 
from hot ring. Black: total spectrum. Red short-dashed: power 
law with F = 1.9. Bottom: Total spectrum divided by power 
law with F = 1.9. 



the surface density of the hot layer is governed by the physics 
of the Coulomb interaction of the hot ions penetrating through 
it, and its temperature by the energy balance between it and 
the underlying disk. With temperature and surface density con- 
strained in this way, the mass flux then depends only on the ra- 
dial drift speed, i.e. the viscosity parameter, a. The actual mass 
flux is low because the temperature of the layer is only about 
80keV. This suggests that the current model is incomplete, and 
we discuss a possible solution in sect.|4] In this paper, however, 
we solve this problem by introducing a parameter, C, by which 
the accretion rate (or equivalently the energy output of each 
component of the flow) is increased. This makes the implicit 
assumption that increasing the accretion rate causes the energy 
output in each component to increase, but the relative contri- 
bution of each component to the overall energy budget and ge- 
ometric configuration of the flow to stay constant. Scaling the 
accretion rate in this way increases the luminosity in all com- 
ponents, which thus causes an increase in temperature for the 
cool disk at a fixed truncation radius. 

With the introduction of the parameters C and f our model 
loses its predictive power, but the goal of this paper is to present 
a model that is plausible rather than precise in its details. In 
the next section we compare our spectra to the best fits from 
observations, and show that for reasonable values of the hot 
layer and inner ring and accretion rate, we can reproduce the 
observed soft excesses using a significantly truncated accretion 
disk. 

To illustrate how our model can be made consistent with 
observed soft excesses, we perform a qualitative comparison 
between the soft ex cesses observed in SWIFT J1753.5- 0127 
(iMilleretal.l l2006a') and GX 339-4 (iMiller et al.1 bOO eb*) and 
our model. For each object we use the estimates for black hole 
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mass, inclination and distance to source presented in those two 
papers and take a = 0.2 as a standard value for the viscosity pa- 
rameter. We then assume a moderate truncation radius and find 
a solution for energy and surface density as was discussed in 
sect.|2] and calculate the spectrum. We change C (the accretion 
rate) and ^ (the ratio between seed photons in the hot layer and 
inner ring) in order to match the luminosity and spectral index 
of the best fit to the observed spectrum, and compare our soft 
excess to the observed one. If necessary we also change the 
amount of interstellar absorption, although for both the cases 
we consider we do not need to change it very much. Given the 
systematic uncertainties in our model, a more statistical com- 
parison to the data is not possible; our goal is instead to demon- 
strate that we are able to reproduce the observed spectra with 
physically reasonably parameters. 

3.1. SWIFT J1 753.5-01 27 

We begin with the spectrum from the source SWIFT J1753.5- 
0127. Miller et al. (2006a) took a 42 ksec XMM-Newton ob- 
servation and estimated an X-ray luminosity (0.5-10 keV) of 
LxILem = 2.6 X m'\dl2,.5 kpc)2(M/10 Mq). They fit the spec- 
trum to a power law with a photon index T - 1.67 (2-10 keV), 
and interstellar absorption of A^h - 2.3 x 10^' cm"^. Fitting the 
spectrum with an absorped power law component alone reveals 
a small soft excess below 2 keV, which they fit to a disk with 
kTin ^ 0.22 keV and R,„ ^ Rs (M/lOMg) (J/8.5kpc)/ cos'^^ /. 

A very truncated disk will be too cool to be observable in 
X-rays, while an untruncated disk will have a higher temper- 
ature than is observed (because of the effects of heating from 
the corona). To compare with the observed spectra we assume a 
moderate truncation radius of 15Rs, which is qualitatively very 
different from an untruncated disk (which will have an inner 
radius between 0.5 - 3Rs depending on the spin of the black 
hole). We make the same assum ptions for mass (M - IOMq) 
and distance {d - 8.5kpc) as in iMiller et alJ (l2006al) and as- 
sume an inclination of ju = 0.5. The simulation then predicts 
an accretion rate of M/MEdd = 1.5 x 10"^ for an efficiency 
of 10%. In order to match the measured flux in the power law 
component, we increase the flux in each component by a factor 
C = 12, to give an accretion rate M/MEdd = 1.7 x 10"^. 

In the hot layer, the predicted surface density profile gives 
the optical depth (r) = 0.87, and the energy balance between 
the cool disk and hot layer (see sect. 12.3b determines a temper- 
ature of kTp^ 75 keV, which gives a photon index of F = 1.87 
(3-10 keV)0 The density of the disk and flux in the hot layer 
also allows us to estimate the ionization parameter, ^ = 50 
erg cm s ', from which we get an albedo and find that the 
Comptonized flux incident on the disk heats it to kT^ =0.11 
keV, which looks like kT^ = 0.19 keV when spectral hardening 
is taken into account. 



'^ The spectral index of a Comptonizing corona will change depend- 
ing on viewing angle, since the mean optical depth will change de- 
pending on whether the disk is face-on or tilted, which will harden the 
spectrum. To calculate the energy in the hot layer and thus its tem- 
perature we take the spectral index for the spectrum integrated from 
0°-180°. 
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Fig. 5. Model spectrum for SWIFT J1753.5-0127, for a trun- 
cated disk with inner radius of 15 Rs with a - 0.2 and 
M - IOMq. Red long-dashed line: Spectrum from the disk 
(modified blackbody plus reflection spectrum). Green dash- 
double dotted line: Spectrum from outer disk. Blue dotted line: 
Comptonized spectrum from hot layer. Orange dash-dotted 
line: Comptonized spectrum from hot ring. Black: total spec- 
trum. Red short-dashed: power law with F = 1.66. 



The predicted surface density in the hot inner ring also al- 
lows us to calculate the optical depth of this layer, which we 
find to be (r) = 0.7. Since we do not have a detailed model 
for the radiative transfer in this component, we take a plausible 
value of kT - 200 ke V, which gives a photon index of F = 1.39 
between 2 and 10 keV. The resulting spectrum is still too soft, 
so we set ( - 0.89, meaning that 89% of the disk and reflection 
photons are used to seed the hot layer, while the rest seed the 
hot inner ring. Figure|5]shows the resulting spectral energy dis- 
tribution for SWIFT J1753. 5-0127. The different components 
shown in the figure are the same as in fig. |4] 

Figure |6] (taken from lMiller et al.l (l2006ah ^ shows the spec- 
trum divided by an absorped power law with an absorption 
column density of A^h = 2.3 x 10^'cm"^. Overlaid we show 
their best fit using the XSPEC "diskbb+pow" model (in red) 
and the soft excess predicted by our model (in green). To ob- 
tain a better fit we change the absorption column density to 
A^H = 2.35 X lO^'cm"^ to render the two excesses effectively 
indistinguishable, even though the temperature and peak fluxes 
of both disks are very different. We discuss the reasons for this 
at the beginning of sect. 14.31 Figure |6] also shows that the devi- 
ation from a simple power law in the range 2-100 keV is less 
than 10%, which is consistent with observation. 



3.2. GX 339-4 

Our second source for comparison is the relatively better con- 
strained X-ray source GX 339-4. This source has an estimated 
mass of 6Mh, distan ce of around 8 kpc and inclination// = 0.9. 
Miller et al.] (l2006bl) took a 280ks XMM-Newton observation of 
this source, which they observed at a flux of about Lx/^Edd - 
0.05(M/10Mo)(D/8kpc). They fit the data with a moderate ab- 



D'Angelo et al.: Soft X-ray components in the hard state of accreting black holes 




10 
Energy [keV] 



100 




1.0 10.0 100.0 1000.0 

Energy [keV] 



Fig. 6. The soft excess in SWIFT J1753.5-0127. The crosses 
show the observed spectrum divided by a simple power law 
fit above 3 keV with a best-fit colum n density of A^h = 
2.3 X lO^'cm^ taken from fig. 2 of IMiUeretal.1 (|2006a). 
The dashed-dotted red line shows the best fit two-component 
XSPEC "diskbb-Hpow" divid ed by the same power law (as re- 
ported in lMiller et al.l (l2006a!) ). The dashed blue line shows the 
ratio between our best fit model (see fig.lD and the same power 



law, with an increase of 0.05x10 
tion. 



,21. 



cm in interstellar absorp- 



sorption column density (A^h = 3.72 X lO^'cm"^) and very hard 
power law (F = 1.44) and find a soft excess below 3 keV 
which they fit to a disk with R,„ = 0.6Rs and kTi„ = 0.38 
keV. Additionally, they observe a broad asymmetric Fe-K line 
with a maximum at 6.9 keV, from which they measure a reflec- 
tion fraction of about 0.2-0.3 and ionization parameter ^ ~ lO-' 
erg cm s"'. Fitting the Fe-K line using relativistic broadening 
suggests an inner radius of 0.7Rs. 

For an assumed radius of the inner edge of the disk Ri„ = 
19.Rs, a good fit for this source is obtained with the following 
set of parameter values: C=77 (or M/MEdd - 0.15), ^=0.63, 
and the temperature of the hot layer and hot ring 75 keV, and 
200 keV respectively. The resulting spectral energy distribu- 
ti on is shown in fig.|7] Figure[8]shows the soft excess observed 
in iMiller et al . (2006b). The luminosity enhancement factor is 



higher than for the SWIFT source because of the higher lumi- 
nosity inferred from the observations. The ionization parameter 
in the reflection region is also higher, with ^ - 550 erg cm s"' . 

Note the structure in the data around 1 keV, which is sug- 
gestive of the structure of the interstellar absorption in this re- 
gion of the spectrum. In fact, we can best reproduce the obser- 
vation for a c olumn A^h = 4.9 x l O^'cm'^, some 30% higher 
than found by lMiller et all (l2006bh . 

As is seen in fig. [8l we also find an Fe-K emission line 
in our fit, with a strength comparable to the observations. The 
combination of parameters that fits the overall spectral shape 
in our model therefore also fits the reflection component of the 
spectrum. Since this component depends on details of the in- 
teraction between the hot layer and the cool disk under it, this 
adds some confidence in this part of the physics of our model. 



Fig. 7. Model spectrum for GX339-4, for a truncated disk with 
an inner radius of l9Rs, a = 0.2 and M = 6Mq. Red long- 
dashed line: Spectrum from the disk (modified blackbody plus 
reflection spectrum). Green dash-double dotted line: Spectrum 
from outer disk. Blue dotted line: Comptonized spectrum from 
hot layer Orange dash-dotted line: Comptonized spectrum 
from hot ring. Black: total spectrum. Red short-dashed: power 
law with r= 1.47. 
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Fig. 8. The soft ex cess in GX 339-4, taken from fig. 3 of 
Miller et al.l(l2006bl) . The crosses show the spectrum divided by 
a simple power law fit to the XMM-Newton and RXTE spectra 
with an absorption column density of A^h = 3.72 x 10^' cm"^. 
The 4-7 keV range, where a broad Fe-K line is clearly visible, 
was omitted from the fit. The dashed blue line shows the ratio 
between our best fit model (see fig. |7]) and the same power law 
fit with interstellar absorption A^h - 4.9 x 10^' cm"^. 

4. Discussion 

In the paper thus far we have produced spectra of an ion- 
bombardment model with a truncated disk for Low-Hard state 
accretion, incorporating a spectrally-hardened blackbody com- 
ponent, physical reflection models from the literature and a 
Monte Carlo Comptonization calculation. 

Detailed fits to observations have shown that the model can 
reproduce the observed soft excess, high-energy spectral index 



10 



D' Angelo et al.: Soft X-ray components in the hard state of accreting black holes 



and the approximate strength of the Fe-K Hne in GX 339-4. We 
have also shown that the structure observed in the soft excess 
in GX 339-4 is consistent with interstellar absorption features, 
suggesting that the absorption column density for this source 
might be underestimated. 

However, we have also found that further work is neces sary 
to develop the flow model outlined in iDeufel et al.l (120021) and 
DS05. Most significantly, a large viscosity in the hot regions 
of the flow (the hot layer and ADAF) is required to increase 
the accretion rate sufficiently to match observations. This may 
suggest the presence of strong ordered magnetic fields in the 
inner regions of the accretion flow, which are also believed to 
be associated with observed jets. 



4.1. Mass flux in the hot layer 

As was discussed in sect. [3] the global accretion rate in the 
present model is limited by the rate at which the hot layer can 
flow over the cool disk. The surface density and temperature 
of the hot layer are in turn narrowly constrained by the physics 
of the Coulomb interaction which allow the layer to form and 
the energy balance between it and the underlying disk. There 
are thus two ways to increase the flow rate in the hot layer: by 
increasing the effective viscosity or providing another mecha- 
nism besides viscous dissipation to transport angular momen- 
tum. 

The interaction between the ion supported ADAF and the 
hot layer provides such a mechanism. The ion supported flow 
is partially suppoted against gravity by gas pressure and rotates 
slower than Keplerian. The mass condensing from the ADAF 
on the hot layer thus acts as a sink of angular momentum, which 
increases the mass flux in the hot layer This effect was not 
included in DS05 and the calculations above. An estimate of its 
importance can be made by evaluating the angular momentum 
exchange a posteriori from the solutions in sect. 3. We find that, 
for the viscosity parameter a = 0.2 assumed for the hot layer, 
the effect increases the mass ffux by a factor 2-3. The effect 
is thus significant, but not sufficient to increase the mass flux 
by the factors indicated by the comparison with observations 
in sect. 3. 

The missing ingredient most likely to lead to the higher 
mass fluxes inferred from the observations may well be 
a strong magnetic field. Strong ordered magnetic fields in 
the inner regions of the flow are implied by the pres- 
ence of jets, especially in the hard X-ray states dis- 
cussed here. A bundle of strong ordered magnetic field 
held together by a disk dBisnovatvi-Kogan & Ruzmaikin 
1 1974 can have field strengths well above those pro- 
duced by magnetorotational turbulence. The angular mo- 
mentum exchange by in t eractio r i of such a bun dle with 



4.2. Spectrum of the hot ring 

The spectrum from the hot ring in our model is hard to pre- 
dict without a more detailed model. The uncertainty lies chiefly 
in the distribution and number of seed photons available for 
cooling. The more photon-starved the hot ring is, the higher 
its temperature will be, and (since the evaporation rate into 
an ADAF scales with T~) the smaller its contribution to the 
overall spectrum. The temperature we assumed for the hot ring 
in practice could be much lower, which would bring it more 
in line with observations of the high energy cutoff observed 
in some spectra (which suggest a maximum temperature of 
about kTg ~ 150 keV). Additionally, the geometric distribu- 
tion of seed photons will change the structure of the Compton 
spectrum. This is because photons that scatter once preferen- 
tially scatter back in the direction they were originally trav- 
elling, and there is a deficit of photons in the first scattering 
hump in the spectrum. This effec t is most p ronounced in the 
plane-parallel case (e.g. Haardt & Maraschil (1991) ), but there 
will also be some anisotropy in the hot ring's spectrum if the 
seed photons are primarily from the disk. We have considered 
only seed photons from the disk, but there may also be photons 
produced from other processes (such as synchrotron emission) 
which would allow the hot ring to cool more efficiently and 
make the effects of anisotropy less pronounced (since the seed 
photons would be travelling through the hot ring in essentially 
random directions). Finally, in this paper we have neglected the 
spectral contribution from the ADAF. Depending on its radia- 
tive efficiency, its contribution could also harden the observed 
high-energy Comptonized spectrum considerably. 



4.3. Comparison with other work 

Our model of a disk truncated at 15-20 Rs and surround- 
ing corona is qualitatively very diffe rent from the untrunc ated 
disk (with Ri„ ~ IR s) models fit by iMiller et al] ( l2006ah and 



the disk dSte hle & Spruit" '2001; 



Igumenshchev et al. 2003 ; 



Naravan et al. 2003 : De Villiers et al. 2005) can be much more 
effective than turbulence parametrized with a viscosity param- 
eter a ~ I. This aspect is beyond the present study, and is a 
promising field for further study. 



Miller et al.l (l2006bl) . and it is natural to ask how the observed 
soft excess can be so small when the radiating area is so much 
larger. The answer lies in several points. The most important 
of these is that the temp erature in our disks is about a factor 2 
smaller than is found by Miller et al.l (l2006al) . so that the flux is 
intrinsically much smaller and (even after the colour-correction 
is applied) most of the flux is cut off by interstellar absorption. 
There is a further reduction from the hot surface layer, which 
upscatters about two-thirds of the photons. Finally, the shape 
of the upscattered photons deviates from a power law at low 
energies, so that measuring the temperature of the soft excess 
depends very sensitively on modelling the Comptonized spec- 
trum correctly. 

The effects of irradiation on the measured truncation ra- 
dius have a lso been studied using a more phenomenological 
approach in iGierhnski et al. J2008b, who re-analyzed the data 
from J1817-330 (JRvkoff et al.l2007 !) to demonstrate that iiTadi- 
ation can increase the measured truncation radius in this source 
(although they assume continual stress at the inner boundary 
of the truncated disk, which we have not done here). In par- 
ticular, they note future plans to test the effects of incomplete 



D' Angelo et al.: Soft X-ray components in the hard state of accreting black holes 



11 



thermalization from incident radiation (cf. sect. 12.4.1b . which 
can further increase the truncation radius. 

Several well-studied sources show some evidence of devi- 



ations from a single power law. iDone et al.l ( l2007h notes that 
spectra from Cyg X- 1 have additional structure in their spectra 
that can be fit with an additional very soft Comptonizing com- 
ponent, while both the source GX 339-4 and Cyg X-1 some- 
times show an excess of very high energy photons compared to 
a fit with a single power law, which suggests a second site for 
Comptonization th at is naturally explairi ed with this model. 

For GX 339-4, iMiller et all (l2()06bl) detected a broad Fe- 
K line, which they fitted with a relativistically broadened 
profile, implying an untruncated disk and a spinning black 
hole. However, the observation of a broadened Fe-K line 
may also be consistent with a truncated disk if the broaden- 
ing instead comes from a n o utflow, as has been suggeste d in 
Done & GierUnskil (|2006|) and lLaurent & TitarchukI IJoOTh . 

The truncated disk picture of hard states in X-ray binaries 
has received support from analyses o f the noise spec t rum o f 
the X-ray variability. In the model by IChurazov et al.l (120011) . 
for example, the characteristic frequency at ~ 0.01 Hz in Cyg 
X-1 coiTesponds to the viscous frequency (~ v,./r) of a geomet- 
rically thin accretion flow at a truncation radius of about ISR^. 
X-ray timing data typically contain several characteristic fre- 
quencies, while theoretical models allow for diff'erent mecha- 



nisms of variability (eg. lGiannios & Spruitj(l2004h '). The model 
presented here does not contain enough physics to make pre- 
dictions about the source and nature of X-ray variability. For 
reference we note that, for a truncation radius at 2QRs, the vis- 
cous frequencies for our model are ~ 0.8 Hz for the hot layer, 
~ 2 Hz for the hot ring, and ~ 50 Hz for the ADAF. 1 Hz is a 
characteristic frequency often observed in Cyg X- 1 . 

The physics of interaction of an ion-supported flow with a 
cool disk is very well defined. The process produces a hot layer 
of tightly constrained thickness and temperature. The physics 
is detailed enough, for example, to be implemented in a nu- 
merical hydrodynamic simulations of the accretion flow. The 
present model falls short of achieving this: the inner hot ring 
where the evaporation into the ADAF takes place, in particular, 
contains parametrizations that would need to be improved with 
a more detailed (2-D) treatment of radiative transfer. Thus the 
model used here has adjustable parameters but compared with 
other hot corona models there is a straightforward path to more 
rigorous calculations. 

Note, however, that the likely presence of a strong ordered 
magnetic field in the inner regions of the accretion (see 4. 1 and 
4.2 above) adds additional physics that is not included either 
in existing hot corona models or the present ion illumination 
model. 

5. Conclusions 

From energetic considerations, the hard spectra observed in the 
low-hard state of LMXBs must be produced by hot (kT^ ~ 100 
keV) matter in the inner regions surrounding the black hole. If 
there is also a much cooler disk present, there will necessar- 
ily be some degree of interaction between the two components, 
and the disk will be somewhat heated by irradiation from the 



hot Com pton izing compo n ent. Th e fits reported in lMiller et al 
(l2006ab and iMiller et all (l2006bl) neglect this interaction by 
fitting the disk and hard component separately. In this paper 
we have shown that incorporating the effects of this interac- 
tion heats the inner regions of a moderately truncated disk so 
that, when coupled with the effects of interstellar absorption, 
the size of the soft excess matches observations. Our work also 
highlights the potential pitfalls of using simple power law or 
analytic Comptonization fits at low energies, which can pro- 
vide significant deviations in the soft X-rays, thus changing the 
shape and intensity of the observed soft excess. 

In the case of GX 339-4, our model predicts an Fe-K com- 
ponent of comparable strength to that observed, although we 
did not do a detailed comparison. However, work by others has 
suggested that part of the broadening in the Fe-K line that was 
observed for GX 339-4 can be attributed to a large outflow, and 
detailed models of Fe-K fluorescence in galactic black holes 
show lines that are much broader than is found in AGN models 
(and which are normally used to fit spectra). 

The model we have envisioned presents several opportuni- 
ties for further improvement, in order to better constrain the in- 
troduced fitting parameters, C and rj. The spectrum from the hot 
ring and ADAF are particularly uncertain, and dependent on a 
more detailed model for the radiative transfer through this re- 
gion, as well as the source and number of seed photons (which 
will set the electron temperature in both regions). The model's 
global accretion rate (which is limited by the rate at which the 
hot layer spills over into the hot ring and then evaporates into 
the ADAF) is also very low, although this can be increased if 
the viscosity in the warm layer can be increased, perhaps as a 
result of accretion through an ordered magnetic field. 

CD'A acknowledges financial support from the National 
Sciences and Engineering Research Council of Canada. 
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